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A  mathematical  formulation  for  the  cathode  of  a  membrane  electrode  assembly  of  a  polymer  electrolyte 
membrane  fuel  cell  is  proposed,  in  which  the  effect  of  unsaturated  vapor  feed  in  the  cathode  is  considered. 
This  mechanistic  model  formulates  the  water  saturation  front  within  the  gas  diffusion  layer  with  an 
explicit  analytical  expression  as  a  function  of  operating  conditions.  The  multi-phase  flows  of  gaseous 
species  and  liquid  water  are  correlated  with  the  established  capillary  pressure  equilibrium  in  the  medium. 
In  addition,  less  than  fully  hydrated  water  contents  in  the  polymer  electrolyte  and  catalyst  layers  are 
considered,  and  are  integrated  with  the  relevant  liquid  and  vapor  transfers  in  the  gas  diffusion  layer.  The 
developed  performance  equations  take  into  account  the  influences  of  all  pertinent  material  properties 
on  cell  performance  using  first  principles.  The  mathematical  approach  is  logical  and  concise  in  terms 
of  revealing  the  underlying  physical  significance  in  comparison  with  many  other  empirical  data  fitting 
models. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is  one  of 
several  possible  candidates  for  the  replacement  of  the  internal  com¬ 
bustion  engine  in  automobiles  to  eliminate  pollutant  emissions. 
The  capabilities  of  PEMFCs  have  been  demonstrated  as  stationary 
residential  and  portable  power  sources  [1].  The  performance  of  a 
PEMFC  depends  on  its  well-designed  membrane  electrode  assem¬ 
bly  (MEA)  to  ensure  optimized  mass  transport  and  electrochemical 
kinetics  under  voltage  or  current  operating  modes.  Since  most  volt¬ 
age  consumption  occurs  in  the  cathode  oxygen  reduction  reaction, 
as  compared  to  the  relatively  facile  hydrogen  oxidation  reaction, 
the  cathode  compartment  of  the  MEA  has  been  the  focus  of  many 
theoretical  and  experimental  studies  [2,3].  Fig.  1  is  a  schematic  of 
the  cathode  part  of  the  MEA,  which  consists  of  a  polymer  elec¬ 
trolyte  membrane,  a  cathode  catalyst  layer,  and  a  cathode  gas 
diffusion  layer.  The  solid  polymer  electrolyte  membrane  plays  the 
role  of  electrode  separator  between  the  anode  and  the  cathode, 
as  well  as  the  proton  transport  vehicle.  Fully  dispersed  platinum 
nano-particles  supported  by  high  surface  carbon  powders  were 
fabricated  in  the  catalyst  layer  to  replace  the  previously  used  plat¬ 
inum  black  so  that  a  greater  active  area  per  unit  volume  is  available 
for  the  electrochemical  reaction  [4].  The  gas  diffusion  layer  (GDL), 
made  of  weaved  carbon  fibers,  is  usually  hydrophobically  treated  to 
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enhance  its  water  repulsion  capability.  The  GDL  facilitates  uniform 
oxygen  mass  transfer  from  the  gas  channel,  and  repulses  water  that 
is  generated  in  the  cathode  catalyst  layer  and  that  is  dragged  from 
the  polymer  membrane  by  electro-osmosis  [  5  ].  The  advancement  of 
cell  performance  is  hinged  on  a  proper  MEA  fabrication  technique 
with  appropriate  physical/chemical  properties. 

Many  mathematical  models  have  been  proposed  as  diagnostic 
tools  to  evaluate  MEA  performance  in  practical  operations.  Two 
conceptually  different  approaches  have  appeared  in  the  litera¬ 
ture  of  PEMFC  1-D  simulations,  that  is,  the  macro-homogeneous 
and  the  heterogeneous  models.  In  the  former,  a  mixture  of  car¬ 
bon  supported  platinum  catalysts  (Pt/C),  polymer  electrolyte,  and 
liquid  water-filled  pores  is  assumed  to  be  in  a  single  homoge¬ 
neous  phase  within  the  catalyst  layer  [4-10].  In  the  latter,  various 
types  of  agglomerates  of  mixed  Pt/C  and  ionomers  with  different 
geometric  dimensions  are  connected  as  a  three-dimensional  net¬ 
work  in  the  catalyst  layer  [11-16].  Among  the  agglomerates,  which 
may  be  covered  with  water,  inter-particle  gaseous  oxygen  diffu¬ 
sion  occurs.  In  addition,  gaseous  oxygen  needs  to  be  dissolved  on 
the  external  surface  and  diffuses  to  the  interior  of  the  agglom¬ 
erates  for  the  electrochemical  reaction  on  the  active  Pt/C  site. 
Both  approaches  have  correlated  well  with  experimental  obser¬ 
vations,  except  that  the  double  Tafel  slope  phenomena  observed 
in  the  cathodic  polarization  is  better  simulated  using  the  agglom¬ 
erate  model  [13].  Computational  fluid  dynamics  techniques  have 
been  developed  for  2-D  and  3-D  simulations  with  considerations 
of  heat  transfer  [17-21],  anisotropic  material  properties  [22-24], 
and  convective  flow  in  the  channel  [25-27]  with  various  degrees  of 
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Nomenclature 

a  effective  platinum  surface  area  per  unit  volume 

(cm-1) 

a0  specific  platinum  surface  area  (cm-1 ) 

A  parameter  in  capillary  pressure  head  expression,  Eq. 

(27),  dimensionless 

Co2  oxygen  concentration  ( mol  cm-3 ) 

c0 2, ref  reference  oxygen  concentration,  defined  as  P/Hq2 

(mol  cm-3) 

C  parameter  in  capillary  pressure  head  expression,  Eq. 

(27),  dimensionless 
dc  catalyst  layer  thickness  (cm) 

dd  gas  diffusion  layer  thickness  (cm) 

dm  membrane  thickness  (cm) 

dw  width  of  the  two-phase  region  within  gas  diffusion 
layer  (cm) 

D  parameter  in  capillary  pressure  head  expression,  Eq. 

(27)  (cm) 

Df_j  binary  diffusion  coefficient  for  i  and  j  species 
(cm2  s-1) 

Dflj  effective  binary  diffusion  coefficient  for  i  and  j 
species  (cm2  s-1) 

effective  diffusivity  of  dissolved  oxygen  in  the  cata¬ 
lyst  layer  (cm2  s-1) 

Dx  water  diffusivity  in  membrane  (cm2  s-1 ) 
e  swelling  expansion  coefficient  of  membrane, 

dimensionless 

fc  parameter  defined  by  Eq.  (6c)  (V-1 ) 

F  Faraday’s  constant  (96,487  C  per  equivalent) 

g  gravitational  acceleration  (980  cm  s-2 ) 

Hq2  Henry’s  constant  for  oxygen  solubility 

(atm  cm3  mol-1) 

io,ref  exchange  current  density  at  the  reference  condition 
(A  cm-2) 

I  cathode  current  density  (A  cm-2 ) 

J0  characteristic  current  density,  defined  by  Eq.  (6f) 

(A  cm-2) 

kp<m  membrane  permeability  (cm2) 

I<  liquid  water  permeability  in  diffuser  at  partial  satu¬ 

rated  condition  (cm2) 

Ki,abs  liquid  water  permeability  in  diffusion  layer  at  fully 
saturation  (cm2) 

mm  membrane  molecular  weight  (g  mol-1 ) 

mw  water  molecular  weight  (g  mol-1 ) 

nd  electro-osmotic  drag  coefficient,  defined  by  Eq.  (21 ), 

dimensionless 

Nj  mole  flux  of  species  i  (mol  cm-2  s-1 ) 

P  total  pressure  (atm) 

A  Pm  pressure  difference  between  two  sides  of  membrane 

in  Eq.  (18)  (gem-1  s-2) 
qw  liquid  water  flux  (cm  s-1) 

R  universal  gas  constant  (8.314 J  mol-1 1<-1 ) 

s  liquid  water  saturation  in  diffusion  layer,  dimen¬ 

sionless 

sc  water  saturation  at  the  catalyst-layer/diffusion- 

layer  interface,  dimensionless 
T  cathode  temperature  (I<) 

vw  molar  volume  of  water  (cm3  mol-1 ) 
vm  molar  volume  of  dry  membrane  (cm3  mol-1 ) 

V  catalyst  layer  potential  (V) 

Vc  cathode  potential  (V) 

V0  open  circuit  potential  (V) 

Xi  mole  fraction  of  species  i,  dimensionless 


x?  mole  fraction  of  species  i  in  channel,  dimensionless 
x?  mole  fraction  of  species  i  at  the  catalyst- 

layer/diffusion-layer  interface,  dimensionless 
x[  mole  fraction  of  species  i  at  the  interface  of  single 

phase  and  two-phase  regions,  dimensionless 
z  coordinate  perpendicular  to  the  face  of  the  gas  dif¬ 

fusion  layer  (cm) 

Greek  letters 

ac  cathodic  transfer  coefficient,  dimensionless 
/3i  parameter  defined  by  Eq.  (6b)  (S  cm-2 ) 

p2  parameter  defined  by  Eq.  (13)  (S  cm-2) 

pm  parameter  defined  by  Eq.  (53)  (S  cm-2) 

s  gas  void  fraction  in  diffusion  layer,  dimensionless 

em  volume  fraction  of  the  ionomer  phase  in  the  catalyst 
layer,  dimensionless 

£0  porosity  of  gas  diffusion  layer  at  zero  water  satura¬ 

tion,  dimensionless 

£w,m  volume  fraction  of  water  in  membrane,  expressed 
by  Eq.  (22),  dimensionless 
0  ionomer  potential  (V) 

(p  parameter  defined  by  Eq.  (6a),  dimensionless 

0  capillary  head  (cm) 

A  parameter  defined  by  Eq.  (48),  dimensionless 

ij  parameter  defined  by  Eq.  (46),  dimensionless 

X  water  content  in  membrane  (mole  of  water  per  mole 

of  sulfonic  groups),  dimensionless 
Xs  saturated  water  content  in  membrane,  dimension¬ 

less 

fiw  viscosity  of  liquid  water  (g  cm-1  s-1 ) 

Pm, dry  density  of  dry  membrane  (g  cm-3 ) 

pw  density  of  liquid  water  (g  cm-3 ) 

Geff  effective  electric  conductivity  of  the  diffusion  layer 
(Scm-1) 

<j m  protonic  conductivity  of  membrane  (S  cm-1 ) 

erff  effective  proton  conductivity  in  the  ionomer  phase 
of  catalyst  layer  (S  cm-1 ) 

rw  water  molar  flux  driven  by  hydraulic  permeation, 
defined  by  Eq.  (30)  (mol  cm-2  s-1 ) 


complexity.  However,  to  facilitate  a  quick  evaluation  of  single  cell 
performance  in  a  typical  laboratory,  it  would  be  better  to  have  a 
simple,  versatile  mechanistic  model  for  the  evaluation  of  MEA  fab¬ 
rication  techniques.  Srinivasan  and  co-workers  [28-30]  proposed 
a  semi-empirical  equation  to  describe  the  polarization  behavior 
of  PEMFCs  that  considers  the  effects  of  electrode  activation  over¬ 
potential  and  membrane  ohmic  resistance.  This  was  the  prototype 
of  many  later  empirical  PEMFC  models  [31-33].  Due  to  the  inabil¬ 
ity  of  Srinivasan’s  equation  to  predict  the  limiting  current  behavior, 
an  additional  correction  factor  was  incorporated  by  Squadrito  et  al. 
[31  ].  Other  modifications  to  account  for  the  mass  transfer  effect  of 
the  GDL  at  high  current  densities  were  made  by  Chu  et  al.  [32]  and 
Xia  and  Chan  [33].  However,  the  proposed  equations  in  those  arti¬ 
cles  lack  sound  physical  grounds  for  a  detailed  description  of  the 
associated  electrochemical  process,  and  the  fitted  parameters  are 
only  applicable  to  a  very  limited  operating  region.  Hsuen  [34-36] 
proposed  a  series  of  mechanistic  models  of  the  PEMFC  that  quantify 
all  significant  potential  losses  that  stem  from  cell  electrochemi¬ 
cal  resistances  and  transport  hindrances.  Specifically,  the  voltage 
drops  are  contributed  to  by  gaseous  mass  transport  in  the  gas  diffu¬ 
sion  layer,  dissolved  oxygen  diffusion  in  the  catalyst  layer,  electron 
transport  of  carbon  fiber  in  the  gas  diffusion  layer,  electron  trans- 
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Schematic  of  the  cathode  part  of  membrane  electrode  assembly  (not  to  scale). 

:  Carbon  fiber 

:  Carbon  supported  Pt  catalyst 

□  :  Electrolyte  phase 

□  :  Liquid  water  phase 

Fig.  1.  Schematic  diagram  of  the  cathode  part  of  membrane  electrode  assembly. 


port  of  the  Pt/C  phase  in  the  catalyst  layer,  the  electrochemical 
reaction  of  oxygen  reduction  on  the  Pt/C  surface,  and  proton  migra¬ 
tion  within  the  polymer  electrolyte  phase.  In  addition,  the  influence 
of  liquid  water  blockage  on  oxygen  diffusion  and  cathode  over¬ 
potential  was  quantified  under  a  fully  humidified  air  feed  [36]. 
The  concise,  analytical  while  mechanistic  formulated  approach  is 
extremely  computationally  efficient.  It  can  be  easily  modified,  if 
desired,  with  a  parameter  estimation  scheme  to  retrieve  useful 
kinetic  properties  of  the  MEA  for  possible  fabrication  improvement. 
The  present  work  is  a  further  extension  of  previous  developments 
to  include  the  case  of  unsaturated  air  feed  in  the  cathode.  It  is 
shown  that  the  position  of  the  water  saturation  front  in  the  gas  dif¬ 
fusion  layer  can  be  explicitly  identified  based  on  the  analytically 
derived  equations.  In  addition,  a  less  than  fully  hydrated  poly¬ 
mer  electrolyte  membrane,  and  the  associated  water  diffusion  and 
electro-osmotic  drag,  are  incorporated  in  the  model.  The  proposed 
mechanistic  model  is  more  versatile  and  realistic  in  terms  of  the 
evaluation  of  kinetic/transport  parameters  for  a  better  membrane 
electrode  assembly  design  when  compared  to  other  empirical  data 
fitting  models  of  limited  applicability. 

2.  Mathematical  model 

The  mathematical  model  for  the  PEMFC  cathodes  considered  in 
the  present  work  is  a  steady-state,  one-dimensional,  and  isother¬ 
mal  model.  The  schematic  diagram  of  the  system  is  illustrated  in 
Fig.  1.  The  system  is  composed  of  a  membrane,  a  cathode  catalyst 
layer,  and  the  adjacent  gas  diffusion  layer.  In  the  catalyst  layer, 
the  macro-homogeneous  assumption  is  applied.  This  assumption 
implies  constant  physical  and  chemical  properties  within  the  layer, 
and  that  catalyst  particles,  in  the  form  of  platinum  clusters  sup¬ 
ported  on  carbon  black,  are  well  mixed  with  proton-conductive 
ionomers.  In  addition,  the  water  content  of  the  ionomer  phase 
within  the  catalyst  layer  is  considered  to  be  invariant  with  posi¬ 
tion.  Since  the  electrochemical  reaction  takes  place  at  the  interface 
between  the  platinum  clusters  and  ionomers,  oxygen  molecules 
have  to  be  dissolved  into  the  ionomer  phase  before  reacting  with 
protons.  Other  assumptions  used  in  formulating  the  model  equa¬ 


tions  are  delineated  as  follows: 


(i)  The  ideal  gas  law  is  applied  for  the  gaseous  species  within  the 
gas  diffusion  layer. 

(ii)  The  total  gas  pressure  within  the  gas  diffusion  layer  is  constant. 

(iii)  The  Stefan-Maxwell  equations  are  used  to  describe  multi- 
component  gas  transport  in  the  gas  diffusion  layer,  and 
the  effective  binary  gas  diffusivities  are  evaluated  using  the 
Bruggeman  expression  [4]. 

(iv)  The  interfacial  surface  area  for  the  liquid  and  vapor  phases  of 
water  in  the  diffusion  layer  is  large  enough  such  that  the  water 
vapor  is  saturated  in  the  presence  of  liquid  water. 

(v)  The  membrane  is  impermeable  to  gaseous  species  and  dis¬ 
solved  oxygen. 

(vi)  The  nitrogen  mole  flux  is  zero  due  to  its  inertness  and  assump¬ 
tion  (v). 

(vii)  The  rate  form  of  oxygen  reduction  follows  a  first-order  expres¬ 
sion  with  respect  to  oxygen  concentration.  In  addition,  the 
reaction  rate  is  described  by  the  Tafel  expression  [37]. 

(viii)  The  phase  equilibrium  of  oxygen  at  the  catalyst- 
layer/diffusion-layer  interface  is  achieved  and  follows  Henry’s 
law  [4]. 

(ix)  Convective  mass  transfer  of  dissolved  oxygen  in  the  catalyst 
layer  is  negligible  due  to  small  values  of  its  mole  fraction  in  the 
mobile  phase,  and  Fields  law  is  applied  to  describe  the  transport 
of  dissolved  oxygen  within  the  catalyst  layer. 

(x)  The  potential  loss  due  to  electron  conduction  in  the  catalyst 
layer  is  considered  to  be  negligible. 

(xi)  The  effects  of  various  flow  field  designs  are  not  included  in  the 
present  model.  That  is,  only  the  mass  transport  normal  to  the 
membrane  electrode  assembly  is  considered. 

2. t.  Model  equations  with  no  liquid  water  in  the  gas  diffusion 

layer  (Case  I) 


Under  the  condition  that  liquid  water  is  not  present  in  the  gas 
diffusion  layer,  the  Stefan-Maxwell  equations  take  the  forms 
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where  Nt  denotes  the  mole  flux  of  species  i,  x*  its  mole  fraction, 
R  the  universal  gas  constant,  P  the  total  pressure,  T  the  cathode 
temperature,  and  an  effective  gas-pair  diffusivity  for  i  and  j 
species  in  the  porous  medium  that  is  evaluated  by 


£)eff  _  £).  .£ 

2-j  “  Ul-Jb 


1.5 

0 


(3) 


in  which  s0  denotes  the  porosity  of  the  diffusion  layer. 

In  the  catalyst  layer,  the  equations  of  mass  conservation  and 
Ohm’s  law  can  be  expressed  in  partially  dimensionless  form  as  [36] 


A~iff  ~  J  (exP v°  -v  +  4>)]}xq2=o 
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The  model  parameters  appearing  in  the  above  expressions  are 
defined  by 


<P  = 


aio,refHo2dc  '  „ 

4 FPDf  ’  1 

°2 


g m  .  f  _  «cF 
dc  '  Jc  RT 


_  Co2  _  Co2 
°2  c02,re/  P/Hq2  ’ 


4FPDt 

°  H02dc 


(6a-f) 


the  catalyst  layer,  the  surfaces  of  the  platinum  clusters  are  effective 
for  the  reaction  only  if  they  are  in  intimate  contact  with  dissolved 
oxygen  and  protons.  In  other  words,  the  platinum  surface,  which  is 
not  covered  by  hydrated  ionomers,  is  considered  to  be  inactive  for 
the  reaction.  Therefore,  as  the  water  content  of  the  ionomer  phase 
within  the  catalyst  layer  is  decreased,  the  amount  of  active  surface 
area  of  the  platinum  clusters  is  expected  to  decrease  as  well.  In  the 
present  work,  a  linear  equation  is  employed  as  a  first  approximation 
to  account  for  such  an  observation,  that  is 


where  ac  is  the  electrode  transfer  coefficient,  a  is  the  effective  plat¬ 
inum  surface  area  per  unit  volume,  io  ref  is  the  exchange  current 
density  at  the  reference  condition,  c02  is  the  dissolved  oxygen  con¬ 
centration  in  the  ionomer  phase,  dc  is  the  catalyst-layer  thickness, 
or^f  is  the  effective  protonic  conductivity  for  the  ionomer  phase, 
is  the  effective  diffusivity  of  dissolved  oxygen  in  the  catalyst 
layer,  F  is  the  Faraday  constant,  P  is  the  cathode  pressure,  V0  is  the 
open-circuit  potential,  V  is  the  catalyst  potential,  0  is  the  ionomer 
potential,  and  Hq2  is  the  Henry’s  constant  for  gaseous  oxygen  and 
its  dissolved  form  in  the  ionomer  phase  at  the  cathode  tempera¬ 
ture.  In  the  above  expressions,  Xo2  represents  the  mole  fraction  of 
gaseous  oxygen  in  the  gas  diffusion  layer,  but  stands  for  the  dimen¬ 
sionless  concentration  of  dissolved  oxygen  in  the  ionomer  phase  of 
the  catalyst  layer.  The  same  notation  is  used  because,  by  such  defi¬ 
nitions,  in  both  regions  the  Xq2  profiles  are  continuous  across  their 
boundary. 

At  the  face  of  the  diffusion  layer  (z  =  dm  +  dc  +  dd),  it  is  assumed 
that  mass  transfer  limitations  are  negligible  and,  thus,  one  has 

*o2  =xbo2  (7) 

xw  =  or  Xn2  =  xft2  (8a, b) 

At  the  diffusion-layer/catalyst-layer  interface  (z  =  dm  +  dc),  this 
requires  that 


Xq2  (catalyst  layer)  =  x02  (gas  diffusion  layer)  (9) 


I0  dx o2 
4  F  dz 

?=° 

dz 


(catalyst  layer)  = 


(10) 

(11) 


The  cathode  potential  (denoted  as  Vc)  is  equivalent  to  the  cata¬ 
lyst  potential  minus  the  ohmic  loss  of  the  diffusion  layer.  Thus,  one 
has 


Vc  =  catalyst) 


I 

P2 


(12) 


and 


h  = 


(13) 


where  aeJ^  is  the  effective  electric  conductivity  of  the  diffusion 
layer  and  dd  its  thickness.  It  is  postulated  that  phase  equilibrium 
is  instantaneously  established  between  the  water  vapor  and  the 
water  content  of  the  ionomer  phase  of  the  catalyst  layer  at  the 
diffusion-layer/catalyst-layer  interface,  which  follows  [7] 

X  =  0.043  +  17.8  (g)  -  39.8(^)2  +36.0(^)3  (14) 


In  addition,  the  value  of  the  effective  protonic  conductivity  for 
the  ionomer  phase,  denoted  as  crff,  is  evaluated  by  [7] 


aef  =  eli5  (0.005147.  -  0.00326)  exp 


) 268  (ro 


(15) 


where  sm  represents  the  volume  fraction  of  the  ionomer  phase  in 
the  catalyst  layer.  For  the  electrochemical  reaction  occurring  within 


(,6) 

in  which  Xs  is  the  value  of  X  at  xw  =  xsw  and  a0  is  the  specific  plat¬ 
inum  surface  area. 

It  is  postulated  that  the  membrane  is  impermeable  to  oxygen; 
thus,  one  has 


dx, 


o2 

dz 


=  0 


(17) 


at  the  membrane/catalyst-layer  interface  (z  =  dm). 

As  stated  above,  the  membrane  is  assumed  to  be  impermeable 
to  gaseous  species.  However,  as  a  pressure  difference  is  applied  to 
the  both  sides  of  the  membrane,  flow  of  liquid  water  is  induced  by 
the  pressure  gradients  within  the  membrane,  which  can  be  quanti¬ 
tatively  described  by  Darcy’s  law.  In  addition,  electro-osmotic  drag 
caused  by  proton  flow  and  back  diffusion  resulting  from  the  gra¬ 
dients  of  water  concentration  also  significantly  contribute  to  the 
overall  water  flux  within  the  membrane.  Accordingly,  the  water 
transport  equation  within  the  membrane  can  be  formulated  by  [36]. 


Nw  +  2  Nq2  —  nd  —  — 


Pm, dry  ^  dX 
mm  1  dz 


6w,mPwkp,m  A  Pm 

mwfiw(  1  +  eX)  dm 


(18) 


where  e  (=0.01 26)  is  the  swelling  expansion  coefficient  of  the  mem¬ 
brane,  according  to  Springer  et  al.  [7],  mm  is  the  molecular  weight 
of  the  membrane,  nd  denotes  the  electro-osmotic  drag  coefficient, 
X  is  the  water  content  in  membrane,  is  the  water  diffusivity, 
kp<m  is  the  membrane  permeability,  dm  is  the  membrane  thickness, 
/  is  the  current  density,  fiw  is  the  water  viscosity,  pw  is  the  water 
density,  mw  is  the  molecular  weight  of  water,  £w?m  is  the  volume 
fraction  of  water  in  membrane,  and  A Pm  is  the  pressure  difference 
between  the  two  sides  of  the  membrane.  The  value  of  pm<drylmm  is 
1  lvm  =  1  / 550  mol  cm-3 ,  where  vm  denotes  the  molar  volume  of  the 
dry  membrane.  The  water  diffusivity  is  evaluated  using  [7] 


DX=DX( 303 1<)  exp 


2416 


(19) 


in  which 


D*(303I<)  =  (-2.86  +  1.897.)  x  10“6,  2  <7.  <  3 


Dx{ 303  K)  =  (6.38  -  1.197.)  x  10“6,  3  <  7  <  4 


D*(303  K)  =  (2.563  -  0.337  +  0.026472  -  0.00067173) 

xlO"6,  4  <7  (20a-c) 


The  electro-osmotic  drag  coefficient  is  calculated  as  [7] 

nd  =  2.5^  (21) 

and  the  volume  fraction  of  water  in  the  membrane  is  estimated 
from  Hsing  and  Peter  [38] 


_  (A  —  1  )vw 

£w,m  —  - - 

Pm  +  Xvw 

where  vw  is  the  molar  volume  of  water. 


(22) 


222 


H.-K.  Hsuen,  K.-M.  Yin  /  Journal  of  Power  Sources  196  (2011)218-227 


At  the  exterior  boundary  of  the  membrane  (z  =  0),  the  boundary  3.  Cathode  potential 
conditions  are  formulated  as 

3.1.  Solutions  for  the  case  without  liquid  water  in  the  gas 
X  =  Xs  and  0  =  0  (23a,b)  diffusion  layer 


2.2.  Model  equations  with  liquid  water  in  the  gas  diffusion  layer 
(Case  II) 


As  liquid  water  appears  in  the  gas  diffusion  layer,  the  diffusion 
layer  can  be  considered  to  consist  of  two  regions,  namely  a  one- 
phase  region  and  a  two-phase  region.  In  the  one-phase  region,  the 
Stefan-Maxwell  equations  derived  above  are  also  applicable.  Since 
the  gradients  of  the  mole  fraction  of  water  vapor  in  the  two-phase 
region  become  vanishingly  small,  the  Stefan-Maxwell  equations 
can  be  further  condensed  to  one  equation,  which  is 


P  dx o2 
RT  ~dz~ 


=  -0-< 


Xo2)JVo2 


X 


1  i  _ _ 

^N2- 02  X°2^w-N2  +  (1  ~Xw  -  x02  )Dq2_w 


(24) 


The  effective  gas-pair  diffusivity  needs  to  be  modified  by  the  water 
saturation,  s,  using  the  Bruggeman  relation 


tf*=Di-jel-5{\-sj1-5  (25) 

In  addition,  the  transport  equations  for  liquid  water  within  the  dif¬ 
fusion  layer  are  also  needed  for  the  overall  water  balance.  Darcy’s 
law  was  employed  to  describe  the  liquid  water  transport  in  the 
diffusion  layer.  Because  it  is  postulated  that  the  variations  of  gas 
pressure  within  the  diffusion  layer  are  negligible,  the  gradients  of 
capillary  pressure  of  the  liquid  phase  become  the  only  driving  forces 
for  the  flow  of  liquid  water.  Capillary  head,  indicated  as  0,  is  the 
driving  force  for  water  flux  qw. 


qw  —  — 


I<{s)Pwg 

Mw 


(_d±\  ds 
y  ds  J  dz 


(26) 


where  I<{s)  is  the  diffusion  layer  permeability  and  g  is  the  gravita¬ 
tional  acceleration. 

diflds  and  I<(s)  are  expressed  as  [39] 

^t  =  -AD[e~A^  +  eA^]  (27) 


The  value  of  current  density,  denoted  as  /,  can  be  related  to  Nq2 
as 

No2=-T  (31) 

The  magnitudes  of  DN2_W  and  Do2_w  only  differ  by  about  7%  due 
to  the  similar  molecular  weights  of  oxygen  and  nitrogen.  With  the 
aid  of  the  approximation  D^_w  =  D^_w  and  Eq.  (31),  Eq.  (1)  can 
be  further  simplified  as 


P  dxw 
RT~dz 


4 FNW  -I  \  N] 

1  Xw 


4fd; 


& 

O  2-W, 


D 


& 

o  2-W 


(32) 


Hsuen  [36]  showed  that  this  approximation  greatly  facilitates  the 
simplification  but  without  introducing  appreciable  errors. 

Eq.  (32)  is  solved  analytically  to  give 


Xw  — 


4  FNW 
4  FNW- 


1~(’ 


-  (  x£  +  xf,  + 


b) 


x  exp 


°2  ^  ^  4HVw^ 
RT 


(z  —  dm  —  dc  —  d(j)~ 


PD\ 


fff 

O  2-W 


( Nw  ~  i) 


(33) 


The  solution  of  Eq.  (2)  can  be  obtained  by  substituting  Eq.  (33)  into 
it,  which  yields 


*o2  =  -*n2  exP 


(z-dm-dc-dd)f 


4  FD\ 


-l  +  N, 


fSf 
o2 — n2 


D\ 


\eff  , 

02  —  w  / 


4FNW- 


+( 


*n.  +*n.  + 


x  exp 


(z  —  dm  —  dc  —  dd) 


o2  _rAN2  T  4 FNw-l) 
RT 


PDeff 

1  U0  2-W 


(34) 


Then,  the  value  of  xw,  denoted  as  x^,  at  the  catalyst-layer/diffusion- 
layer  interface  (z  =  dm  +  dc)  is  readily  calculated  as 


*w  — 


4  FNy, 


4FNW- 


J-  (xo2+xn2  +  4FN’w  _  , ) 


exp 


-ddRT 

PDeff 

1  U0  2-W 


(35) 


I<(s)  =  Khabs(s  +  0m )  (28) 

where  Ki  abs  is  the  absolute  diffusion  layer  permeability.  The  overall 
water  balance  within  the  diffusion  layer  requires  [36] 

[(2  +  4nd)JVo2  +1 Vw  +  rw]  (^)  +  (|)  =  0 

(29) 

and  rw  denotes 

^  _  Pwkp,m(^s  —  l)rwAPm  (30) 

mwflw  (vm  +  ^sVw)  (1  +  ^s)dm 


For  fixed  values  of  /  and  Nw,  Eqs.  ( 1 8)-(23a)  form  an  initial  value 
problem  with  one  state  variable  X.  In  the  present  work,  the  integra¬ 
tion  of  Eq.  ( 1 8)  from  the  membrane  exterior  boundary  (z  =  0)  to  the 
catalyst-layer/membrane  interface  (z  =  dm)  is  carried  out  using  an 
explicit  Euler  method.  After  the  integration,  /  remains  unchanged 
while  Nw  is  adjusted,  and  the  above  integration  procedures  are 
repeated  until  the  value  of  X  at  z  =  dm  matches  the  value  estimated 
using  Eqs.  (14)  and  (35)  within  the  desired  tolerance.  Since  only 
one  unknown,  Nw,  is  to  be  solved,  a  bisection  method  is  employed 
for  iterations.  As  the  value  of  Nw  is  determined  for  a  given  /,  the 
values  of  x02  and  xN2  at  the  diffusion-layer/catalyst-layer  interface 
(z  =  dm  +  dc),  indicated  as  xc0  and  x^2 ,  are  then  calculated  as 


At  the  boundary  between  these  two  regions,  this  requires  that 
xw  =  xsw  as  well  as  the  continuation  of  Xo2  and  Xn2-  The  governing 
equations  for  the  catalyst  layer  remain  the  same  except  that  the 
water  content  of  the  ionomer  phase  reaches  its  saturated  value. 
In  addition,  the  water  transport  equation  for  the  membrane  is 
no  longer  needed  since  the  profile  of  the  water  content  within  it 
becomes  uniform. 
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*N2  =*N2eXP 


-ddRT 


-I 


Nn 


4  FDe^ 

KHru02-N2 


fff 

O 2-W  , 


D\ 


Inserting  Eqs.  (25)  and  (43)  into  Eq.  (44)  yields 


(37) 


dx07 


3.2.  Solutions  for  the  case  with  liquid  water  in  the  gas  diffusion 
layer 

Under  the  condition  that  liquid  water  appears  in  the  gas  diffu¬ 
sion  layer,  the  width  of  the  two-phase  region,  denoted  as  dw,  can 
be  directly  calculated  using  Eqs.  (18),  (31 )  and  (33)  by  assigning  the 
mole  fraction  of  water  vapor  at  the  boundary  of  the  one-phase  and 
the  two-phase  regions  to  its  saturated  value,  which  yields 


(1  -Xsw-Xo2) 

where 


=  -rj 


(s  +  0.01)  [e-^s-c)  +  e^s-c)]  ds 


(1-s)1 


(45) 


q  = 


Kl,  absRTADgp \ 


1 


Pmw/xwSQ5  [2  +  An^  +  *w/(l  —  *w)  —  4rwF//j  _^n2-o2 

_ _ 

( 1  —  )^02  -vi 


+ 


(46) 


4  FPDeft 
Hri  uo 2-W 


[(4nd  +  l)/-4Frw](l  -xbw)+I 


a  _rj _ kj2-w _ 1  1  l v  '  *  r  11  ■'wjv1  w /  1  J  l  c  r 

w~  d  RT [(4nd+l )/-4Frw]  }  [(4nd  +  l)/-4Frw] (1  -xsw)  +  /  f  xo2  =  *o2  “*n2  IexP(^)-  1] 


After  integrating  Eq.  (45)  over  the  two-phase  region,  one  eventually 
arrives  at 


(38) 

The  mole  fraction  of  Xq2  and  Xn2  at  the  boundary  of  these  two 
regions,  denoted  as  and  xj^ ,  are  readily  determined  using  Eqs. 
(34)  and  (33),  and  the  constraint  of  1  =  x^2  +  *02  +  xjy,  which  give 


*02  =  -*n2  exP  \  -(dd  -dw)^p 


-I  +(4nd  +  2)/-4Frw 


Deff 

.  °2-n2 


yff 
o  2-w 
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O2  '  n2  ^  (4nd  +  l)/-4Fru 


x  exp  -(dd-dw) 


FT[(4nd  +  l)/-4Frw 


4FPD02-w 


(39) 


\  -(dd-dw)RT 

-I  t  (4nd  +  2)7  -  4Frw 

1 

4  FP 

Deff  +  Deff 

lU02-N2  uO  2-w  J 

J 

(2 

>i2  =  *n2  exP' 


In  the  two-phase  region,  the  Stefan-Maxwell  equations  can  be 
rearranged  from  Eq.  (1)  to  give 


Nw  —  Nq2 


=  Nn 


Dn2-w*o2  +  Dq2-w  (l  -  X5W  -  Xo2 ) 


(41) 


With  the  approximation  Def  =  D/J  Eq.  (41 )  is  further  simpli- 


tied  to 

n-Mt ^5) 


N  2-w 


fff 


0  2-w’ 


(42) 


Inserting  Eq.  (42)  into  Eq.  (29),  the  overall  water  conservation 
within  the  diffusion  layer  results 


(2  +  4nd)  + 


£*N2 -w*w 


Ni 


02 


TTw  — 


f^N2-w^02  +  ^02-w(l  ~  ~  *02) 

Khobsp2wgAD(s  +  0.01 )  +  e^)]  /  ds  \ 


/xwmw 


E  -°  (43) 


Using  the  same  approximation  D^_w  =  D^_w,  and  inserting  Eq. 
(31 ),  Eq.  (24)  gives  the  form 


o2  -  ^o2  “  *n2 

in  which 

;  0  (s  +  0.01 )  [e-^s-c)  +  e^s-c)] 


-[ 


ds 


(47) 


(48) 


(1-s) 

where  sc  represents  the  water  saturation  at  the  diffusion- 
layer/catalyst-layer  interface. 

3.3.  Performance  equations 

When  an  unsaturated  air  stream  is  employed  as  a  cathode  feed 
for  a  PEMFC,  the  amount  of  water  transported  from  the  exterior 
membrane  boundary  and  generated  by  the  oxygen  reduction  is  usu¬ 
ally  not  enough  to  fully  hydrate  the  membrane  and  the  ionomer 
phase  in  the  catalyst  layer  at  low  current  densities.  The  same  sit¬ 
uation  might  be  encountered  for  a  saturated  air  feed  if  a  positive 
pressure  difference  is  imposed  between  the  two  sides  of  the  mem¬ 
brane.  Under  such  a  condition,  liquid  water  is  not  present  in  the  gas 
diffusion  layer,  and  consequently  the  equations  for  case  I  should  be 
employed  for  the  calculations.  As  the  current  density  is  increased 
to  a  value  beyond  which  the  ionomer  phase  in  the  catalyst  layer 
can  be  completely  hydrated,  liquid  water  starts  to  appear  in  the 
diffusion  layer.  For  such  an  instance,  the  equations  for  case  II  are 
used  in  the  computations. 

After  the  water  content  profile  within  the  membrane  is  cal¬ 
culated,  the  protonic  conductance  of  the  membrane  is  readily 
computed  using  Eq.  (15).  For  the  case  that  liquid  water  emerges 
in  the  diffusion  layer,  the  membrane  is  fully  hydrated;  thus,  one 
has 

1 


am  =  (0.005139AS  -  0.00326)  exp 


,268 (303 


t) 


(49) 


If  liquid  water  is  absent  in  the  diffusion  layer,  the  value  of  am  is 
estimated  as  an  averaged  value  using 

rdm  1 

; 


CTm  — 


=  —  r 

dm  Jo 


(0.0051 39A.  -  0.00326)  exp 


1268 


V303 


i) 


dz 

(50) 


Once  the  values  of  xcQ ^  x^,  and  om  become  available,  the  cathode 
potential  can  be  determined  using  the  performance  equations  [36] : 

v  =  Vr  __L _ L _ L_1 

C  °  2/1,  h  Pm  fc 


x  In 


I  +  P\  //c[exp(-J/c//J, )  -  1  ] 


<p[I0xc02  -  P\  /fc  +  exp(-//c/2/J,  )(/J,  /fc  -  I0xc02  + 1/2)] 

(51) 
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l2[2xCo2I0fc  -  Pt  +  Pi  exp(-2x^2/o/c/^i )]  1 

Mxco2)2lZ[xco2I0fC  ~Pi+Pi  exp (-xc02I0fc/fa )]  J 

for  /  >  2/0Xq2  (52) 

in  which 

Pm  =  ~7—  (53) 

Q-m 

The  cathode  potential  expressions  from  Eqs.  (51)  and  (52)  can 
be  viewed  as  the  open-circuit  value  minus  the  individual  potential 
losses  stemming  from  the  limitations  of  reactant  transport,  proton 
migration,  electron  conduction,  and  the  electrochemical  reaction. 
They  were  derived  by  assuming  that  the  oxygen  and  ionomer 
potential  profiles  within  the  catalyst  layer  are  parabolic  polyno¬ 
mials  (Eq.  (51))  or  piecewise  parabolic  polynomials  (Eq.  (52)).  The 
switch  between  these  two  equations  should  be  made  as  the  oxy¬ 
gen  depletion,  computed  using  the  prescribed  profiles,  occurs  at  the 
catalyst-layer/membrane  interface.  The  procedures  of  their  deriva¬ 
tion  and  the  expressions  for  the  individual  potential  losses  will  not 
be  presented  here;  readers  interested  in  the  details  of  the  equation 
formulation  are  advised  to  consult  the  listed  literature  [36]. 

4.  Results  and  discussion 

The  performance  equations  developed  in  the  preceding  sections 
were  used  to  construct  the  discharges  curves  of  PEMFC  cathodes 
so  as  to  investigate  the  effects  of  the  model  parameters  on  the 
cathode  performance.  Since  the  emphasis  of  the  present  work  is 
on  exploring  the  influence  of  water  saturation  within  the  diffu¬ 
sion  layer  on  the  cathode  performance,  only  the  analyses  of  certain 
model  parameters  whose  changes  may  give  rise  to  particular  char¬ 
acteristics  of  liquid  water  transport  within  the  diffusion  layer  will 
be  included  here.  The  parameter  values  listed  in  Table  1  are  used  in 
the  calculations  unless  otherwise  specified.  It  should  be  noted  that 
the  formulas  and  the  parameter  values  for  the  capillary  pressure 
of  the  gas  diffuser  were  taken  directly  from  the  work  of  Natara- 
jan  and  Nguyen  [39]  without  any  modifications.  Different  forms 
of  capillary  pressure  based  on  other  considerations  have  been  also 
reported  in  the  literature  [40,41  ].  As  shown  in  the  procedures  of  the 
model  formulation,  the  present  approach  is  quite  general  and  is  not 
restricted  to  any  particular  forms  of  capillary  pressure  of  the  gas  dif¬ 
fusion  layer.  Consequently,  other  expressions  of  capillary  pressure 
can  also  be  incorporated  into  the  present  performance  equations 
without  appreciable  difficulties. 

Fig.  2  demonstrates  the  discharge  curves  for  cathode  feeds  with 
different  relative  humidities.  Higher  vapor  content  indeed  favors 
an  earlier  saturation  in  the  membrane  and  catalyst  layer  so  that 
better  proton  conductivity  is  expected;  as  a  result,  high  vapor 
saturation  outperforms  in  the  low  current  density  region.  For  a 
better  illustration  of  the  hydration  effect,  the  emergence  of  liq¬ 
uid  water  in  the  gas  diffusion  layer  is  marked  in  the  figure.  This 
point  indicates  the  initial  occurrence  of  a  fully  hydrated  membrane 
and  catalyst  layer,  before  which  the  polymer  membrane  and  the 
adjacent  catalyst  layer  are  only  partially  hydrated.  The  associated 
current  density  is  determined  by  letting  dw  =  0  in  Eq.  (38),  and 
the  oxygen  mole  fraction  between  one-phase  and  two-phase,  x £  , 
can  be  calculated  directly  from  Eq.  (39).  This  oxygen  concentration 
is  equivalent  to  the  interfacial  concentration  of  CL/GDL  x^.  The 
corresponding  cathode  voltage  Vc  is  then  determined  either  from 
Eq.  (51)  or  (52)  based  on  the  criterion  imposed  on  the  calculated 
current  density.  Note  that  for  RH  =  1,  full  hydration  commences 
from  the  open  circuit  potential  for  a  steady-state  operation.  Before 
membrane  saturation,  ohmic  resistance  varies  and  decreases  with 
discharged  current  until  full  hydration  is  arrived.  As  the  membrane 
is  saturated,  liquid  water  flooding  begins  at  the  interface  of  the 


Fig.  2.  Polarization  curves  at  various  cathode  feed  relative  humidities.  dm  =  50  (Jim, 
anode  relative  humidity  =  1 ,  cathode  gas  diffusion  layer  porosity  =  0.5,  diffusion  layer 
liquid  water  permeability  =  3  x  10_10cm2,  pc  =  pa  =  5  atm,  temperature  353  K.  Solid 
symbols  coordinates  (0.3452,  0.7109),  (0.6760,  0.6397)  are  conditions  of  emerging 
flooding  in  cathode  gas  diffusion  layer  for  cathode  feed  relative  humidities  RH  =  0.5 
and  0,  respectively. 

catalyst  layer/gas  diffusion  layer.  Liquid  water  in  the  gas  diffusion 
layer  obstructs  the  gaseous  mass  transfer  of  oxygen  from  the  chan¬ 
nel.  Drier  cathode  feed  results  in  a  higher  limiting  current  density 
due  to  a  favorable  gas  diffusion  of  oxygen  through  the  gas  diffusion 
layer. 

The  water  saturation  profiles  corresponding  to  various  cur¬ 
rent  densities  at  the  specified  inlet  RH  =  0.5  are  depicted  in  Fig.  3. 


dimemsionless  distance 

Fig.  3.  Water  saturation  profiles  in  the  gas  diffusion  layer  at  different  discharge 
current  densities  for  a  specified  cathode  feed  RH  =  0.5.  Other  conditions  are  the  same 
as  that  in  Fig.  2. 


Table  1 

Values  of  model  parameters. 
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Membrane  thickness,  dm  (cm) 

Catalyst  layer  thickness,  dc  (cm) 

Gas  diffusion  layer  thickness,  dd  (cm) 

Effective  electric  conductivity  in  diffusion-layer  and  catalyst-layer,  creJf  (S  cm-1 ) 

Effective  ionomer  conductivity  in  catalyst  layer  in  fully  hydration,  crjf  (S  cm-1 ) 

Effective  diffusivity  of  dissolved  oxygen  in  the  catalyst  layer,  (cm2  s-1 ) 

Gas-pair  diffusivity,  Do2-n2  (cm2  s-1 ) 

Gas-pair  diffusivity,  Dn2_w  (cm2  s-1 ) 

Gas-pair  diffusivity,  Do2_w  (cm2  s-1 ) 

Swelling  expansion  coefficient  of  membrane,  e 

Porosity  of  gas  diffusion  layer  at  zero  water  saturation,  e0 

Volume  fraction  of  the  ionomer  phase  in  the  catalyst  layer,  em 

Product  of  platinum  surface  area  and  reference  exchange  current  density,  aoi’o.re/  (A cm-3) 
Gas  pressure,  P  (atm) 

Saturated  vapor  mole  fraction  at  353  K,  xsw  (atm) 

Bulk  oxygen  to  nitrogen  molar  ratio,  /x^ 

Temperature,  T  (K) 

Henry’s  constant,  Ho2  (atm  cm3  mol-1 ) 

Open  circuit  potential,  V0  (V) 

Cathodic  transfer  coefficient,  ac 

Hydraulic  permeability  of  membrane,  kPtm  (cm2) 

Molar  volume  of  dry  membrane,  vm  (cm3  mol-1 ) 

Molar  volume  of  liquid  water,  vw  (cm3  mol-1 ) 

Water  viscosity,  /jlw  (gem-1  s-1) 

Density  of  liquid  water,  pw  (g  cm-3 ) 

Liquid  water  permeability  in  diffusion  layer  at  fully  hydration,  K[Abs  (cm2) 

Constants  in  capillary  heat  expression  A,  C,  and  D  (cm) 

Fully  hydrated  water  content  in  membrane  (mole  of  water  per  mole  of  sulfonic  groups),  ks 


0.005 
0.0005 
0.03 
1  [36] 

3.1157  x  10-2 
1  x  10-5  [36] 

5.58  x  10-2 
7.74  x  10-2 
7.40  x  10-2 
0.0126 
0.5 
0.5 

1  x  10-2  [36] 

5  [36] 

0.0934 
21/79 
353  [36] 

2.04  xlO5  [4] 

1.2  [4] 

1.0  [36] 

1.58xl0-14  [4] 

550 

18 

3.565  xlO-3 
0.972 

3  x  10-10  [36] 

3.7,  0.494  and  0.0173  [39] 
14.043  [7] 


As  shown,  the  water  saturated  domain  in  the  gas  diffusion  layer 
expands  along  with  the  increased  current  density  until  the  lim¬ 
iting  condition  is  arrived.  Since  the  assumption  of  instantaneous 
phase  equilibrium  is  imposed  in  the  model  formulation,  it  is  also 
desirable  to  discuss  the  characteristics  of  the  inter-phase  mass 
transfer  of  water  molecules  within  the  two-phase  domain.  Remem¬ 
ber  that  for  a  fixed  current  density,  Nq2  remains  constant  over  the 
entire  diffusion  layer,  as  exhibited  by  Eq.  (31).  When  the  original 
Stefan-Maxwell  equations  are  applied,  under  the  presence  of  liq¬ 
uid  water,  the  values  of  Nq2  and  Nw  are  related  by  a  function  of 
xsw,  Xo2,  and  gas-pair  diffusivities  as  revealed  by  Eq.  (41 ).  This  indi¬ 
cates  that  small  amounts  of  liquid  water  evaporate  in  the  two-phase 
region  since  the  value  of  Xq2  decreases  monotonically  towards 
the  catalyst-layer/diffusion-layer  interface.  Nevertheless,  the  sit¬ 
uation  changes  as  the  approximation  D^_w  =  D^_w  is  employed. 
Under  such  a  condition,  Nw  remains  unchanged  in  the  two-phase 
region,  as  implied  by  Eq.  (42).  Consequently,  evaporation  of  liquid 
water  only  occurs  at  the  water  saturation  front.  After  liquid  water 
vaporizes,  a  major  portion  of  it  continues  to  travel  along  the  same 
direction  to  the  gas  channel,  while  the  rest  is  dragged  by  the  flow 
of  oxygen  to  the  catalyst-layer/diffusion-layer  interface.  It  should 
be  stated  that  in  the  work  of  Natarajan  and  Nguyen  [39],  a  small 
number  (0.01)  was  added  to  the  water  saturation  in  Eq.  (28)  so  as 
to  alleviate  the  infinitely  large  slopes  at  the  water  saturation  front 
and  to  circumvent  the  numerical  difficulties  encountered  in  their 
computations.  However,  such  an  addition  is  not  necessary  in  the 
present  study  since  calculations  with  Eqs.  (47)  and  (48)  are  still 
feasible  without  the  presence  of  this  small  number. 

The  individual  potential  losses  caused  by  the  limitations  of 
oxygen  electrochemical  reduction,  oxygen  mass  transport,  pro¬ 
ton  migration,  electron  conduction,  and  membrane  ohmic  drop 
for  cathode  feeds  of  relative  humidities  0  and  1  are  presented  in 
Fig.  4.  In  general,  the  activation  overpotential  associated  with  oxy¬ 
gen  reduction  increases  sharply  as  the  current  density  is  driven 
from  the  open  circuit  voltage.  The  rate  of  increase  is  alleviated  as  the 
current  density  is  further  raised.  As  a  comparison,  a  higher  activa¬ 
tion  overpotential  is  predicted  when  dry  air  is  fed  instead  of  air  with 
saturated  vapor.  The  membrane  and  the  catalyst  layer  are  more  eas¬ 


ily  hydrated  if  a  more  humidified  cathode  feed  is  employed,  and 
the  activation  overpotential  is  reduced  due  to  a  greater  catalytic 
surface  area  being  available  for  the  electrochemical  reaction.  High 
inlet  humidity  also  reduces  the  ohmic  loss  in  the  membrane  and 
protonic  ohmic  loss  in  the  catalyst  layer  due  to  greater  hydration 


Fig.  4.  Polarization  curves  of  individual  contributive  terms  calculated  using  per¬ 
formance  equations  at  cathode  feed  relative  humidities  0  and  1.  (♦  0)  Activation 
overpotential  for  oxygen  reduction  reaction,  (tv)  diffusion  overpotential  of  cata¬ 
lyst  layer,  (a  A)  diffusion  overpotential  of  gas  diffusion  layer,  (•  O)  ohmic  potential 
loss  of  ionomer  phase  in  catalyst  layer,  (■□)  ohmic  potential  loss  of  gas  diffusion 

layer,  and  ('/At'aT  )  ohmic  potential  loss  in  the  membrane.  Open  symbol  for  RH  =  1, 
filled  symbol  for  RH  =  0.  dm  =  50  |jim,  anode  relative  humidity  =  1 ,  cathode  gas  diffu¬ 
sion  layer  porosity  =  0.5,  diffusion  layer  liquid  water  permeability  =  3  x  10-10cm2, 
Pc  =  Pa  =  5  atm,  temperature  353  K. 
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Fig.  5.  Effect  of  membrane  thickness  and  pressure  difference  between  cathode  and 
anode.  Cathode  pressure  Pc  =  5  atm.  Anode  relative  humidity  =  1,  cathode  feed  rel¬ 
ative  humidity  =  0.5,  cathode  gas  diffusion  layer  porosity  =  0.5,  temperature  353  K. 
Solid  symbols  indicate  the  current  densities  of  initial  flooding  in  gas  diffusion  layer. 
(1.445,  0.5223),  (0.6201,  0.5754)  are  for  AP=2atm,  membrane  thickness  50  and 
200 |jim,  respectively.  (0.3452,  0.7109),  (0.3452,  0.6694)  are  for  AP  =  0atm,  mem¬ 
brane  thickness  50  and  200  p,m,  respectively. 

in  the  ionomer  phase,  as  indicated  in  Fig.  4.  The  major  part  of  the 
overall  ohmic  loss  of  the  cathode  is  caused  by  the  proton  migration 
in  the  ionomer  phase  of  the  catalyst  layer  and  electron  conduction 
in  the  gas  diffusion  layer.  The  potential  loss  in  the  gas  diffusion  layer 
is  linearly  proportional  to  the  discharged  current  density.  Oxygen 
is  depleted  in  the  cathode  catalyst  layer  at  high  polarization,  and 
the  associated  diffusion  overpotential  increases  dramatically  as  the 
limiting  current  density  is  approached. 

The  effect  of  membrane  thickness  and  the  pressure  difference 
between  the  cathode  and  anode,  A P,  is  demonstrated  in  Fig.  5. 
An  overall  better  performance  is  achieved  when  using  a  thinner 
membrane  as  a  result  of  less  ohmic  resistance,  better  water  diffu¬ 
sion  rate  in  the  membrane,  and  lower  water  saturation  in  the  gas 
diffusion  layer.  The  pressure  effect  differs  at  the  low  and  high  dis¬ 
charged  regions.  Fligher  values  of  A P  prohibit  the  hydration  of  the 
membrane  in  the  early  stage,  but  prevent  water  flooding  in  the  gas 
diffusion  layer  at  high  overpotential.  Fig.  5  also  shows  that  the  back 
pressure  is  more  influential  if  a  thinner  membrane  is  used  because 
a  greater  pressure  gradient  is  imposed.  For  instance,  for  a  specified 
AP  =  2  atm,  a  50-p.m  membrane  becomes  fully  hydrated  at  as  high 
as  1.445  A  cm-2,  in  contrast  to  that  of  0.6201  A  cm-2  if  a  200-[xm 
membrane  is  used.  The  corresponding  water  saturation  profiles  at 
the  limiting  current  condition  are  plotted  in  Fig.  6.  A  higher  back 
pressure  gradient  prohibits  the  convective  water  flow  and  induces 
less  water  flooding  in  the  cathode  gas  diffusion  layer. 

Fig.  7  demonstrates  the  effect  of  void  fraction  in  the  gas  diffusion 
layer  on  the  polarization  curve  at  RFI  =  0.5.  It  shows  that  the  mem¬ 
brane  and  the  catalyst  layer  are  saturated  at  a  higher  overpotential 
if  a  gas  diffusion  layer  of  greater  porosity  is  used.  There  is  also  a 
slight  performance  enhancement  in  the  low  current  density  region 
for  the  diffusion  layer  with  a  low  porosity  due  to  better  hydration  of 
the  membrane  and  the  catalyst  layer.  As  liquid  water  floods  exten¬ 
sively  within  the  diffusion  layer,  an  MEA  fabricated  with  larger 
gas  pores  allows  for  easier  evaporation  of  liquid  water  and  facil¬ 


Fig.  6.  Water  saturation  profiles  at  limiting  current  densities  at  various  membrane 
thickness  and  pressure  difference  between  cathode  and  anode.  Cathode  pressure 
Pc  =  5  atm.  Key  is  the  same  as  that  in  Fig.  5. 

itates  better  water  vapor  diffusion;  as  a  consequence,  better  cell 
performance  results. 

Fig.  8  evaluates  the  effect  of  liquid  water  permeability  in  the  gas 
diffusion  layer,  J<j  abs,  on  the  cell  performance.  Large  values  of  I<iabs 
induce  a  facile  convective  mass  transport  of  water,  as  suggested  in 
the  third  term  on  the  left  hand  side  of  Eq.  (43);  in  consequence, 
lower  water  saturation  is  expected  in  the  porous  medium.  Oxygen 


I  /A  cm'2 


Fig.  7.  Effect  of  gas  diffusion  layer  void  fraction  on  the  cell  performance. 
dm  =  50  |jim,  anode  relative  humidity  =  1,  diffusion  layer  liquid  water  permeabil¬ 
ity  =  3  x  1 0-10  cm2 ,Pc  =  Pa  =  5  atm,  temperature  353  K.  Coordinates  (0.2470, 0.7398), 
(0.3452,  0.7109),  and  (0.4538,  0.6849)  are  the  conditions  of  emergent  flooding  in 
gas  diffusion  layer  for  porosities  0.4,  0.5,  and  0.6. 
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I  /A  cnrv2 


Fig.  8.  Effect  of  diffusion  layer  liquid  water  permeability  on  the  cell  performance. 
dm  =  50  |jim,  anode  relative  humidity  =  1,  cathode  gas  diffusion  layer  porosity  =  0.5, 
Pc  =  Pa  =  5  atm,  temperature  353  K.  (0.3452,  0.7109)  is  condition  of  emergence  of 
flooding  in  gas  diffusion  layer. 

transfer  is  more  favorable  when  less  gas  void  is  blocked  by  the  liq¬ 
uid  water  in  the  diffusion  layer,  as  is  most  prominent  in  the  high 
current  density  region.  Note  that  before  the  emergence  of  liquid 
water  at  the  membrane/catalyst-layer  interface,  cell  performance 
is  unaffected  by  the  parameter  I<i<abs. 

5.  Conclusions 

Mechanistic  performance  equations  of  unsaturated  cathode 
feed  are  derived  for  the  cathode  part  of  the  MEA  of  a  general  PEMFC. 
Individual  potential  losses  within  the  MEA  can  be  assessed  explic¬ 
itly  from  first  principles  without  the  physical  ambiguity  that  is 
encountered  in  empirical  data  fitting  models.  In  the  model,  water 
transport  in  the  membrane  is  described  with  electro-osmotic  drag, 
diffusion,  and  pressure  induced  convection  mechanisms  in  either 
partially  or  fully  hydrated  situations.  Detailed  oxygen  and  vapor 
concentrations  in  the  diffusion  layer,  water  contents  in  the  ionomer 
phases  within  the  membrane  and  catalyst  layers,  and  the  liquid 
water  saturation  profile  in  the  diffusion  layer  are  all  well  integrated 
in  the  analytical  formulation.  The  interaction  of  Stefan-Maxwell 
type  ternary  gaseous  diffusions  and  liquid  water  transport  are  cor¬ 
related  by  the  capillary  equilibrium  correlation  within  the  porous 
medium.  The  model  simulates  PEMFC  cathode  performance  from 
dry  to  fully  humidified  cathode  feeds.  Oxygen  concentration  and 
liquid  water  saturation  distributions  in  the  gas  diffusion  layer  are 


explicitly  determined  for  the  specified  current  density  and  humidi¬ 
fication.  Computation  of  the  performance  equation  is  efficient,  and 
is  applicable  to  an  extensive  operating  range,  which  allows  for  a 
better  MEA  design  with  desired  physical/chemical  properties. 
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